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Design  of  Energetic  ionic  Liquids  (PREPRINT) 

Jerry  A.  Boat:.  Air  Forte  Research  Laboratory .  Space  and  Missile  Propulsion  Division. 

Edwards  AFB.  CA:  Hut  Li.  Department  of  Chemistry.  University  ofNebraska-Llncoln ;  Mark  S. 
Gordon,  Department  of  Chemistry.  Iona  Stale  University. 

An  essential  need  of  the  IJS  Air  Force  is  the  discovery,  development,  and  fielding  of  new, 
energetic  materials  for  advanced  chemical  propulsion  in  space  and  missile  applications.  Some  of 
the  key  factors  driving  the  requirement  for  new  chemical  propellants  include,  (a)  improved 
performance  in  terms  of  increased  specific  impulse  and  density,  (b)  reduced  sensitivity  to 
external  stimuli  such  as  impact,  friction,  shock,  and  electrostatic  discharge,  and  (cl  mitigation  of 
environmental  and  toxicological  hazards  (and  the  resulting  costs)  associated  with  currently  used 
propellants. 

A  class  of  compounds  which  can  potentially  meet  tlicse  requirements  is  known  as  ionic  Liquids 
(ILs).  w  hich  are  chemical  salts  with  unusually  low  melting  points.  The  physical  and  chemical 
properties  of  ILs  render  them  useful  for  many  purposes,  most  notably  as  environmentally  benign 
(“green")  solvents/reaetion  media  but  also  as  catalysts,  electrolytes.  etc.[  1 )  From  a  Department 
of  Defense  (DoD)  perspective.  ILs  ate  being  explored  as  new  propellants  and  explosivcs.[2]  The 
Air  Force,  in  particular,  is  interested  in  ILs  as  potential  replacements  for  currently  used 
monoptopellants  such  as  hydrazine  which  is  carcinogenic,  highly  toxic,  and  has  relatively 
modest  performance  characteristics.  In  contrast,  many  ILs  have  superior  densities  and  specific 
impulses  as  well  as  significantly  reduced  sensitivity  and  toxicity  characteristics.  Furthermore, 
their  properties  can  he  carefully  tuned  via  the  choice  of  the  component  ions. 

The  overall  objective  of  the  Design  of  Energetic  Ionic  Liquids  challenge  project  is  to  address 
several  key  technical  issues  and  challenges  associated  with  the  characterization,  design,  and 
development  of  ILs  as  new  monopropellants.  Among  these,  for  example,  are  a  fundamental 
understanding  of  the  (instability  of  ILs.  the  intrinsic  nature  of  the  short-  and  long-range  structure 
and  interactions  between  the  component  ions|2e-f|.  and  identification  of  live  key  steps  in  the 
initial  stages  of  decomposition  and  eombustion|2a-c].  The  research  described  in  this  article  is 
focused  on  characterization  of  the  structures  and  stabilities  of  ion  pair  clusters  and  prediction  of 
their  interaction  energies  in  the  gas  phase. 

Our  computational  approach  utilizes  quantum  chemical  metltods  for  prediction  of  ion  pair 
structures  and  interaction  energies.  In  particular,  geometry  optimizations  were  performed  using 
second-order  perturbation  theory(3]  (MP2.  also  known  as  MBPT<2))  with  the  aug-ce-pvdz  basis 
set[4|.  denoted  as  MP2/aug-cc-pvdz.  Relative  energies  were  refined  using  a  systematic  scries  of 
single-point  energy  calculations  at  the  MP2  and  coupled  cluster  (e.g..  singles  and  doubles  with  a 
perturbative  estimate  of  triples.  CCSD(T)[S1)  levels  of  theory.  Specifically.  MP2'cc-pvdz. 
MPl'aug-cc-pvdz.  and  CCSD(T)'ce-pvdz.  energy  calculations  were  combined  to  obtain  estimated 
CCSDtTy aug-ce-pvdz  relative  energies.  All  computations  were  performed  using  the  CiAMESS 
quantum  chemistry  code.|6J 

MP2  and  coupled  cluster  (CC)  calculations  in  GAMF.SS  utilize  a  library  of  communications 
routines  known  as  the  Distributed  Data  Interface  <DDI).[7]  a  high-level  communications  layer 
operating  between  CiAMESS  and  the  underlying  message-passing  protocols  (Shared  Memory 
(SHMEM).  Message  Passing  Interface  (MPI).  Low-level  Application  Programming  Interface 
(L.API).  or  sockets  within  a  Transmission  Control  Protocol; Internet  Protocol  (TCP/IP)  stack.)  In 
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IBM  systems  KRAKEN  and  BABBAGE.  DDI  uses  MPI  for  intranode  communications  and  lire 
LAPI  protocol  for  messages  between  nodes.  These  types  of  calculations  have  significant  memory 
requirements  and  therefore  are  well  suited  for  execution  on  systems  w  ith  large  amounts  of 
memory  per  node,  such  as  BABBAGE.  Coupled  cluster  calculations  are  especially  memory 
intensive  and.  as  implemented  in  GAMESS  with  DDI.  utilize  a  three-fold  hierarchy  of  memory. 
First,  a  modest  amount  of  Replicated  Data  (RD)  is  exclusively  assigned  to  each  core.  Similarly,  a 
block  of  Node-specific  Data  (ND)  is  reserved  on  each  node  and  is  shared  by  all  live  cores  on  that 
node.  The  remaining  memory  on  each  node  is  collectively  shared  by  all  cores  as  a  large,  single 
pool  of  Distributed  Data  ( DD>.  Therefore,  the  required  Memory  ( Mcc I  per  node  for  CC 
calculations  is  M, ,  =  P*(RD|  *  (NDl  +  (DD)W,  where  "P"  and  "N"  are  (l>c  number  of  cores  per 
node  and  the  total  number  of  nodes,  respectively,  used  in  the  computation. 

The  values  of  RD.  ND.  and  DD  are  determined  by  the  specifics  of  the  calculation,  whereas 
suitable  values  of  P  and  N  are  dictated  by  the  hardw  are;  specifically,  the  amount  of  accessible 
physical  memory  per  node.  If  necessary.  P  can  be  chosen  to  be  smaller  than  the  number  of 
available  cores  per  node  Pm«  in  order  to  reduce  the  amount  of  required  memory  per  node.  Table 
I  summarizes  the  memory  requirements  for  CCSD.T)  calculations  using  a  series  of  increasingly 
large  basis  sets.  Only  the  smallest  calculation  (CCSDiTVcc-pvdzl  could  be  performed  w  ithin  the 
constraints  of  tire  hardware  (Pmm  and  Masa<)  and  the  challenge  queue  limits  (Nnui  and  Tm«.  sec 
Table  2.)  In  principle,  the  CCSDiT|/6-3 1  l«-*G(d.p)[8]  and  CCSIX Tl  'aug-ce-pvd/  calculations 
could  be  run  on  the  pair  of  "bigmem"  nodes,  but  the  estimated  required  wall  time  of  live  former, 
on  the  order  of  100  days,  is  prohibitively  long.  Conversely,  this  calculation  would  be  within  the 
realm  of  practicality  if  — 100  bigmem  nodes  w  ere  available. 

One  of  the  specific  ion  combinations  considered  in  this  w  ork  is  the  1 ,2,4-triazolium  cation 
([C.N.H.l’i  paired  w  ith  the  dinitramide  anion  ||N(NO-)>J-).  Of  the  numerous  structures  found 
for  the  two  paiis  of  1 ,2,4-triazolium  and  dinitramide  ions,  or  the  pairs  of  corresponding  neutral 
1 ,2,4-triazolc  and  dinitramine  molecules,  the  most  stable  MP2'aug-cc-pvdz  optimized  geometries 
are  shown  in  Figure  I .  In  the  ionic  structure,  each  1 ,2,4-triazolium  forms  two  hydrogen  bonds, 
via  the  hydrogens  on  the  N  atoms,  to  tlie  O  atoms  of  the  dinitramide  ions.  Interestingly,  this 
structure  exhibits  parallel  stacking  of  the  two  cationic  1 ,2,4-triazolium  rings.  The  interplane 
distance  is  -3.2  A.  with  a  parallel  displacement  of  - 1 .4  A.  The  corresponding  neutral  tetramcr 
shows  a  similar  parallel  stacking  arrangement  of  the  triazole  rings. 

Furthermore,  it  is  of  interest  to  determine  the  cluster  size  at  w  hieh  the  ion  pair  structures  become 
more  stable  than  the  corresponding  neutral  pair  structures.  A  previous  study  predicted  that  ion 
pair  dimers  are  typically  higher  in  energy  than  neutral  pair  dimcrs.[2c|  Including  zero  point 
vibrational  energy  (ZPV’E)  corrections,  the  ionic  tetranter  in  Figure  I  is  1.2  kilocalorie  /  mole 
(kcal/mol)  lower  than  that  of  the  neutral  one.  The  MP2  method  tends  to  picdict  higher  energies 
for  ionic  species  vs.  neutral  speeies.(2e)  so  more  accurate  CCSD(T>'aug-cc-pVDZ  energy 
calculations  of  these  two  tetranter  structures  were  desired.  However,  since  the  computational  cost 
of  CCSD<T)i'aug-cc-pVDZ  is  prohibitive,  tltese  energies  were  approximated  from  tlte  MP2/aug- 
cc-pVDZ  energies  by  estimating  the  electron  correlation  energy  differences  using  three 
independent  methods:  (1)  the  differences  between  the  MPl'cc-pVDZ  and  CCSDTy'cc-pVDZ 
energies  of  the  tetramers.  (2)  the  differences  between  the  MPI'aug-ec-pVDZ  and  CCSDfTyaug- 
cc-pVDZ  energies  of  the  twelve  pairs  of  dimers  in  these  two  tetramers.  and  (3)  the  differences 
betw  een  the  MP2aug-cc-pVDZ  and  CCSD(Tlaug-ce-pVDZ  energies  of  the  eight  monomers  in 
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these  two  tetramen.  Using  these  tluee  methods,  and  including  ZPVE  corrections,  the  estimated 
CCS D(T >aug-cc-p V DZ  energy  of  tlie  ionic  tetramer  is  low  er  than  that  of  the  neutral  tetramer  by 
5.7, 7.3.  and  7.7  kealmol,  respectively. 

In  conclusion.  quantum  chemical  calculations  suggest  that  cation-cation  parallel  stacking 
structures  can  exist  in  scry  small  ionic  clusters  such  as  two  1,2,4-triazolium  cations  and  two 
dinitiamide  anions.  Furthermore,  for  two  pairs  of  1 ,2.4-triazolium  and  dinitiamide.  ionic 
structures  are  more  stable  than  the  corresponding  neutral  structures.  Finally,  it  should  be  noted 
that  lower  theoretical  mctltods  which  do  not  include  the  efi'ccts  of  electron  correlation,  such  as 
Hanree-Fock.  do  not  predict  a  parallel  stacking  geometry  of  the  rings.  Therefore,  it  is  essential  to 
utilise  correlated  methods  such  as  MP2  ami  CCSD(Tl  in  order  to  obtain  proper  descriptions  of 
the  structures  and  interaction  energies  of  these  ion  clusters.  The  structural  motifs  and  interaction 
patterns  found  in  this  study  provide  new  under  standing  of  ionic  materials  with  aromatic  rings. 
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Tabic  I.  Memory  requirements  for  CCS1XT)  single  point  energy  calculations. 


Basis  set 

Iff  of  AOs) 

KD 

(MW.’ core) 

ND 

(MW/node) 

DD 

(MW) 

i 

N 

MCC 

(MW:  node) 

cc-pvdz 

s 

1175 

4950 

64 

1-381 

(376) 

m 

6-31 1— t-G(d.p) 

22 

3298 

16.000 

64 

3.900b 

(580) 

64 

3.570b 

19 

2 

ll.474c 

aug-cc-pvdz 

26 

3.875 

19.150 

64 

4.200b 

<622) 

1 

2 

13.476c 

aug-cc-pvtz 

330 

18.493 

146.000 

64 

21. 105b 

<1268) 

a 

2 

91.823d 

aug-cc-pvqz 

2495 

60.370 

807.000 

64 

75.475b 

(2228) 

n 

2 

466J65d 

a  “MW"  denotes  megawords  (106  64-bit  words.) 

b  Exceeds  amount  of  usable  physical  memory  on  each  standard  node  (see  Table  2.) 

c  Fits  within  usable  physical  memory  on  each  bigmem  node,  but  execution  time  is  prohibitively 
long. 

d  Exceeds  amount  of  usable  physical  memory  on  each  bigmem  node  (see  Table  2.) 
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Tabic  2.  Challenge  and  bigmem  queue  characteristics  on  tabbage. 


Uueuc 

Coics  pet  node 

(Pmax) 

Maximum  #  o!  nodes 
(Nmax) 

Maximum 
memory  node  Mmax 

(MW/node) 

Maximum 
wall  time 

Tmax 

<hnuis) 

eluillenge 

16 

(vl 

-3,500 

48 

16 

2 

-15.500 

48* 

a  -MW"  denotes  megawords (106  64-bit  words.) 

*  48  hour  limit  obtained  via  special  request.  The  default  wall  time  limit  of  the  bigmem  queue  is 
1 2  houni. 


Figure  I .  M  Pl'a  ug  -  ce-p  V  DZ  optimized  structures  of  two  pairs  of  1 ,2,4-triazolium  1 1 .2.4- 
triazole)  and  diiutruimde  idiiuUainine)  molecules.  II  is  white.  C  is  gray,  O  is  red.  N  is  blue 
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